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We study the dynamics of strings by means of a distribution function /(A,B,x, t) defined on 
a 9+ ID piiase space, wliere A and B are tiie correlation vectors of riglit- and left-moving waves. 
We derive a transport equation (an analogous to Boltzmann transport equation for particles) that 
governs the evolution of long strings with Nambu-Goto dynamics as well as reconnections taken 
into account. We also derive a system of coupled transport equations (an analogous to BBGKY 
hierarchy for particles) which can simultaneously describe long strings /(A, B, x, t) as well as simple 
loops /(A,B,x, f) made out of four correlation vectors. The formalism can be used to study non- 
linear dynamics of fundamental strings, D-brane strings or field theory strings. For example, the 
complicated semi-scaling behavior of cosmic strings translates into a simple solution of the transport 
system at small energy densities. 

PACS numbers: 98.80.Cq 11.27.+d 11.25.Yb ll.25.-w 



I. INTRODUCTION 



The semiclassical evolution of long fundamental strings I ^ R ^ Is (where R is the size of compact dimensions 
and Ig is string length) is described fairly accurately by Nambu-Goto action until the strings start to intersect in the 
target space. Such intersection can lead to reconnections with non vanishing probability p > which significantly 
complicates the dynamics. The problem proved to be relevant not only for fundamental strings or D-brane strings, 
but also for field theory strings on various gravitational background. In fact most of numerical simulations of strings 
were designed to study cosmic strings with p 1 TUT] . The numerical studies revealed a very interesting semi-scaling 
phenomena: scaling of large scales and non-scaling of small scales [SHlOj which one can understand reasonably well 
even analytically [TltiTS] . 

In addition to numerical simulations a couple of attempts where made to describe the dynamics of strings by means 
of a transport equation on a probability distribution of only string lengths [121 HZ] above Hagerdon temperature or 
on a probability distributions of string intersections with a given two dimensional surface [H] . Unfortunately, not all 
of the transport phenomena can be described in terms of these probability distributions. In this paper we describe 
the statistics of strings with a probability distribution of the left- and right-moving correlation vectors (representing 
diamond-shaped section of the world-sheet) which are directly related to observable quantities. The phase space 
remains reasonably small (only nine dimensional) but all of the relevant transport phenomena (such as Nambu-Goto 
evolution, transverse and longitudinal reconnections, production of loops and background gravitational effects) can 
be included in the dynamics. 

The paper is organized as follows. In the next section we define correlation vectors and construct the corresponding 



phase space. In Sec. Ill the transport equation is derived for long strings with Nambu-Goto evolution and transverse 
collisions taken into account. The dynamics of a system with long strings and simple loops is analyzed in Sec. |IV| 
The main results are summarized in the conclusion. 

II. CORRELATION VECTORS 

The classical dynamics of strings in the limit of zero coupling gs = is described by Nambu-Goto action: 

^NG = -i^ J ^ -dei{hab)dadt (1) 

where hab — 9fiuX'^ax'^b- ^'-'^ convenience we also assume that 27ra' — 1 and Ig = ^^/l^- The corresponding equation of 
motion is 
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After gauge fixing (/iqi = 0, /iqo + hu = 0, and t = x^) on a fiat background we get a system of equation for 

three- vectors 

X - x" = 0, x^ + x'2 = 0, x' • X = 0, (3) 

whose solution can be described in terms of left- and right-moving waves 

^ a(a-t) + b(a + 

with condition |a'| = |b'| = 1. 

In general the exact evolution of a given infinitesimal segment at co is known only if all of the functions sl{(t) and 
b(a) for all strings are known. In a kinetic theory such evolution should be describable by a system of equations 
analogous to the BBGKY hierarchy for interacting particles. Truncating the system at a giving order is equivalent 
to neglecting higher order correlations which is often a starting point in the kinetic theory. At the leading order the 
truncation of BBGKY hierarchy for particles leads to the molecular chaos assumption such that interacting particles 
are uncorrelated. Following the same logic for strings we define an ensemble of possible continuations of a(cr) and 
h{a) away from ao described by correlation vectors 



A(cro) = -( f daa'{ao + a)) and B(c7o) = ( / dah\ao + a)) (5) 

\J-L/2 / \J-L/2 I 



for sufficiently large L^^=|A|,L^i?=|B| (Note a minus sign in the definition of A.) The exact value of L 
is unimportant and one can safely send L — ^ oo for the correlation vectors defined on infinite strings. Describing the 
ensemble of continuations with correlation vectors and ignoring higher statistical momenta is an assumption analogous 
to the molecular chaos assumption for particles. We will refer to it is as the strings chaos assumption. 

The string chaos assumption is equivalent to approximating left- and right-moving waves by piece-wise linear 
functions such that neighboring segments are completely uncorrelated. Then the world-sheet can be viewed as a 
collection of diamond shaped regions glued together to form a curved surface and an arbitrary point with correlation 
vectors A and B would be located somewhere on the diamond (A, B). In the limit where the string coupling is small 
(but non-zero) the Nambu-Goto action gives a good description of the dynamics of long strings Z ^ i? ^ up until 
the point when strings start to pass through each other in the target space. A possible outcome of such intersection 
includes reconnections of right- and left-moving waves with non-zero probability whose exact value would depend 
on compacification details. The probability is typically small for fundamental or D-brane strings p <C 1 but can be 
quite large for field theory strings p « 1. In either case the system with a large number of strings becomes quickly 
intractable, but one might hope to understand the dynamics by means of a kinetic theory. If one applies the strings 
chaos assumption to the entire network and approximates all of the strings by piccewise linear functions then the 
entire world-sheet of all strings including all reconnections would be approximated by a collection of diamonds. Then 
the main problem is to describe the statistical evolution of the system with a transport equation (analogous to the 
Botzmann transport equation for particles). 



III. TRANSPORT EQUATION 

Let us define the energy density (or invariant length density) of strings /(A, B, x, t) with correlation vectors A and 
B at spatial position x and time t. Then the evolution of this distribution from time t till time t + 5t\s given by 

m «+m 



/(A, B, X + v( A, B)«, ( + «) = /( A, B, X, «) + ^ ^ j « + (^^ J « + (^^ J ^ it + ^ j ^ « (6) 

where the expected velocity of an infinitesimal segment due to the string chaos assumption is v(A, B) = | ^A -|- B^ 
and A = ^ , B = ^ are the unit vectors. The corresponding transport equation is 

where the gravitational term is expected to be given by some difi'erential operator 



dt), 



-H/ (8) 
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which depends on the background geometry (see Appendix |A|) . The dependence of /, v and H on A, B and x will be 
suppressed throughout the paper for brevity of notations. The three "collision" terms on the right hand side of ([t]) 
represent: Nambu-Goto evolution, transverse (i.e. intersection of nearby strings) and longitudinal (i.e. production of 
simple loops) reconnections. The exact role played by these terms will become more transparent once we derive them 
below. 

To simultaneously describe the populations of long strings as well as small loops it will be convenient to define two 
distribution functions / and / such that / = / + /, however at very large energy densities it is expected that most 
of the energy is in long strings / ~ /. In this limit we can neglect the production of small loops « in (7 1: 

■ v-l^+H/^l^l +1^1 , (9) 



dt 9x 




NG 




and study the remaining terms (Nambu-Goto and transverse) one by one. 

As the time passes the neighboring world-sheet diamonds exchange correlation vectors due to classical dynamics 
described by the Nambu-Goto action. Since the world-sheet area of a given diamond D (see Fig. [T]) with correlation 




FIG. 1: Diamond D with all of its nearest neighbors W, E, N, S, NW, NE, SW and SE. 



vectors A and B is AB/2, the probability of a random infinitesimal segment on D to exit D through NW boundary 
in unit time is ^-^-^ — l/A and the probability to exit D through NE boundary is 1/B. Similarly a given 
infinitesimal segment on ether SW or SN diamonds can enter D through SW and SN boundaries respectively. In 
the transport equation such evolution would be describe by two destruction and two creation terms representing 
transitions of infinitesimal segments through NW, NE, SE and SW boundaries respectively, i.e. 

(¥) - - + /.A-.B-^'(^ + /dA-.B-^'(^ (10) 

A B J AA'Na J BB'Nb ^ ' 



where 



iV../f?^/(A,B) «,,/:^/,A,B). (II) 
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The transverse collisions describe interactions of diamonds which could be arbitrary far away on the world-sheet but 
nevertheless intersect in the target space. Such intersection might produce sudden deformations of the world-sheet 
with non-zero probability p> 0. See Fig.M The area spanned by a diamond (A, B) is ^\AaB\, where A= (1, 




FIG. 2: Diamond undergoing a transverse collisions at point O. 



B = (1, are the null four-vectors. Therefore the probability for an infinitesimal segment to cross MN or OP 
{PQ or NO) per unit time is (x 1/A {(x 1/B) and the probability to collide with an infinitesimal segment of another 
diamond (A',B') is proportional to a four- volume given by \{A A B) A {A' A B')\. All of these transitions (across 
MN, NO, OP and PQ boundaries) lead to modification of the correlation vectors to be described by a collision 

term 

completely uncorrelated suggesting the following transverse term 



(f)^ m the transport equation. Under the strmgs chaos assumption the colliding diamonds are taken to be 



5/' 
dt 

P 



-p I dA'dB'^^^-^liA A B) A {A' A B')\f{A,B)f{A',B') 



0<a,l3<l 



dadA'dp<m'^-^±^ (\{A A B') A [A' A i3)|/(A', f )/(-, B') + |(i A B) A (i' A f )/(A', B') 

4 \ pa a p 



where the integration over a and /3 describes an unknown location of the intersection on participating diamonds . 
The first and second lines in (12 1 represent respectively the destruction and creation of infinitesimal segments due to 
transverse collisions. 



It follows that the transport equation ([t]) together with Nambu-Goto ( 10 ) and transverse ( 12 ) terms takes the 
following simple form: 



dadA'dpdB' 

{A + B 



1 



1 



AA'Na BB'Nb 

AabaA'ab' 



(12) 



/(A',B)/(A,B')- 



(13) 



/(A',f)/(-,B') 

p a 



/(^,^)-/(A,B))/(A',B') 



This equation can already be used to study networks of long strings at large energy densities in Minkowski space 
(where H = 0), but it is also possible to derive H for more complicated background geometries. In the Appendix [a| 
we provide a sample derivation of the operator H for Friedniann backgrounds which is most relevant for cosmology. 
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IV. HIERARCHY OF EQUATIONS 



Until now we have only concentrated on long strings, but it is straightforward to generalize our discussion to 
include simple loops. One can show that the simplest loops must contain at least two right-moving and two left- 
moving correlation vectors: Ai, A2, Bi and B2 where Ai -I- A2 = Bi -|- B2 and Ai + A2 = Bi + B2 is the total 
invariant length of the loop. Such loops are describe by four diamonds (Ai,Bi), (Ai,B2), (A2,Bi) and (A2,B2) 
that would overlap in the target space if the left- and right-moving waves were peciewise linear made out of {Bi, B2} 
and {Ai, A2} respectively. One might argue that this corresponds to an infinite number of self-intersections, but as 
we will see the longitudinal collisions of long strings produce non-self-intersecting loops which are only approximated 
with four correlation vectors under the string chaos assumption. 

More generally one might also want to include populations of more complicated loops consisting of 5, 6 and more 
correlation vectors with additional distribution functions and additional transport equations. A complete system of 
such equations would be analogous to the BBGKY hierarchy for particles, but we shall truncate the hierarchy at only 
two equations describing long strings / and the simplest loops / made out of 4 correlation vectors: 



Of 




at 


ax 


Of 




dt 


ax 



+ fl/-(S)„, + (§l), + (S), 



(14) 



This should be sufficient for the purpose of understanding the leading as well as next-to-leading order phenomena at 
small energy densities. 



To describe the Nambu-Goto evolution of correlation vectors on long strings /(A,B) we can still use (10), but the 
dynamics of correlation vectors on simple loops /(A, B) should be slightly modified, i.e 



dt 



/(A,B) /(A,B) 



NG 



A 



B 



dA'dB' 



1 



1 



AA'Na BB'Nb 



/(A', B)/(A, B')5(A+A'-B-B') W^'-B-B') (15) 



where 



iVA(A,B) 
^b(A,B) 



dA'dB' ^ 

-^1 Jy-^ j-f jO(A+A'-B-B')"(A+A'-B-S') 

dA'dB' ^ 

^1 /(A i-t> jC>(A+A'-B-B')"(A-|-A'-S-B')- 



(16) 



Note that in contrast to the worldsheet of long strings described by (10), the world-sheet of loops is not generated 
with completely random diamonds but only with diamonds which can be on the same simple loops enforced by two 
delta functions (5(a+A'-b-B') and 5(a+A' -b-B')- 

The transverse collisions between diamonds lead to sudden deformation of the worldsheet regardless of whether the 
participating diamonds belong to long strings or simple loops. Almost surely the aftermath of such collisions would 
be strings that are not describable with only four correlation vectors, resulting in a flow of energies from / to /. Thus, 
the transverse term derived above ( 12 ) can be easily generalized to simultaneously describe populations of simple 
loops 



as well a long strings 



Ah B A A' A 13' 



/(A',B')/(A,B) 



(17) 



df 
dt 



p / dA'dB 



>{A + B) 



AaISaA' AB' 



/(A',B')/(A,B) 



(18) 



p / dadA'dpdB 



,{A + B) 



AaB aA' AB' 



/(A',|)/(^,B') + /(^,f)/(A',B') 

pa a p 



Note that the second line in ( 18 ) represents creation of infinitesimal segments and involves an integral of the energy 
density of all strings / = / + /• By substituting ( |lo} [isj [is} [l7| in (14) we obtain a system of transport equations. 
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d_ 

at 



Ox 



Ox 



H 



A 



^AA'Na 



H 



A^B^A'^B' (/(a', f )/(^,B') + (/(^, f ) - /(a,b)) /(a',b')^ ^ „^ , 

/ = /dA'dB' + ^) /(A',B)/(A,B')W-B-B')5(A+A'-B-B') 



^)/(A',B)/(A,B': 

B 

/3 



BB'Nb 
{A+B) ■ 



A A B A A' A B' 



/(A',B')/(A,B) 



where only two longitudinal terms at") and (^§|) remain to be calculated. 

In contrast to the transverse collisions which involve interactions of only two diamonds (or four correlation vectors) , 
the longitudinal collisions must involve at least three nearby diamonds described with at least five correlation vectors. 
This is due to the fact that any pair of diamonds is allowed to undergo a longitudinal collision only if they are not 
nearest neighbors (e.g. E and NW diamonds on Fig. Tjcan collide, but E and D cannot). In particular the collisions 
between next-to-nearest diamonds might produce a loop which is describable with only four correlation vectors in 
terms of the distribution function /. For example, the simplest longitudinal collisions takes place when a pair of left- 
moving correlation vectors collides with a triplet of right- moving correlation vectors or vise versa (see Appendix |b] 
for details). Altogether the longitudinal production of loops is given by (Bl), or 



dt 



Ji Y[ da,dAidAdBj/(A„B,) J (5(5-^^^^^(„^A, 



,1=1,2,3 



)'5(a2-l)%) 



A. 



'^(A-A2)'^(B-/3iB2 



B2 



(20) 



(A o B) 



where (A o B) denotes the same as a previous term but with all of A's (i.e. A, Ai, A2 and A3) interchanged with all 
of B's (i.e. B, Bl, B2 and B3) . In addition to the longitudinal production of simple loops, the longitudinal collisions 
lead to a merger of correlation vectors on long strings resulting in growth of correlations described by ( B3|B4 ). Both 
integrals can be combine into the following longitudinal term 



(21) 



^(Bi-B2)'^(a2-l)'^(fe) 



'5(A-(l~Qi)Ai-(l-a3)A3) X! '^(B-(l-ft)BO i I - Pi + {1 - Oi) 



i=l,3 



1=1,3 



- '5(A-A,)'5(B-B.) ( 1 + ]^ ) ^ '5(A-A2)'5(B-B2) 



(A^B) 



This concludes our derivation of the transport phenomena for strings. The system of transport equations ( 19 ) 



together with longitudinal (20 21 ) terms determine the leading as well as next to leading order dynamics of strings. 



Although the system might seem a bit complex, one can solve it fairly easily either numerically or perturbatively in 
various limits. In fact one such limit was already studied in |14 | ll5 j . where the semi-scaling dynamics of cosmic strings 
was first analyzed. One can show that the transport equation for long strings ( 19 1 reduces to an equation similar to 



Eq. (17) of Ref. [T3] in the limit where the longitudinal production of the smallest loops dominates the non-linear 
dynamics on long strings. 



V. CONCLUSION 



The main objective of the paper was to derive transport equations which could describe a very complicated non- 
linear evolution of strings. We have shown that at large energy densities the transport phenomena can be describe 
with a single equation ([t]), analogous to the Boltzmann transport equation for particles. The key ingredient in the 
construction was the distribution function /(A,B,x, i) defined on a 9-l-lD phase space, which is sufficiently large to 
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allow all of the relevant transport phenomena such as Nambu-Goto evolution, transverse collisions and gravitational 
effects. 

At small energy densities the analysis was further complicated by a necessity to simultaneously describe populations 
of long strings and small loops. This was accomplished by deriving a system of transport equations which is analogous 
to BBGKY hierarchy for particles. By truncating the hierarchy we have concentrated on a pair of coupled equations 
( 19 1 which is sufficient to describe leading as well as next-to-leading order transport phenomena (including longitudinal 
collisions) for long strings /(A,B,x, i) as well as simple loops /(A,B,x, i). 

The next step should be a careful analysis of the solutions of the transport equations, but it is already clear that 
the formalism developed in the paper should be useful in describing complicated non-linear systems with strings (e.g. 
fundamental strings, D-brane strings or field theory strings). For example the analysis of cosmic strings, which is 
usually tackled numerically due to non-trivial interaction between small and large scales, should be greatly simplified 
with the help of derived equations. 
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Appendix A: Friedmann background 

The simple decomposition of the string solution into left- and right-moving waves Q is only valid in Minkowski 
space, but on other gravitational backgrounds the situation is more complex. As an example, consider the string 
motion x((t, t) on a Friedmann background 

ds^ = a{tf{dt^ - d^^) (Al) 

described by Hubble rate = - in conformal coordinates. One can show that when two straight segments of opposite 
moving waves pass through each other the Hubble friction would produce interactions resulting in correlations of 
previously uncorrelated segments. Thus, our main task is to estimate the effect of such interactions on the correlation 
vectors A and B in comoving coordinates. 
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From the equation of motion ^ the tangent vectors satisfy 



-a'^i/(b'-(a'.b')a') 
|b' = H(a'-(b'.a')b') 
and the comoving string length density evolves as 

(See Ref. [H] for details.) Then the correlation vectors defined in ^ must obey 

A = -h/J ' (b'(CT) - (a'((T) • b'((7))a'(cr))\ = -CiHB - C2H(1 - 2v2)A 



and 



B 



-U / y"^^' (a'(a) - (b'(c7) • a'(a))b'(a))\ = -C^HK - C^H{1 - 2v^)B 



(A2) 
(A3) 

(A4) 
(A5) 

(A6) 



where the factors Ci, C2, C3 and C4 can be estimated as follows. The first term — CiffB represents an action of the 
correlation vector B which is according to the integral of ( A5 ) must takes place even if the segments A and B do not 
overlap. In other words B must be acting on the correlation vector A for it entire duration ^ A even if segment B 
acts on a given infinitesimal segment only during time interval ^ B. To take this effect into account we set Ci = A/B 
or similarly C3 — B/A. In contrast to the first term, the second term — C2-ff(l — 2v^)A does not vanish only when 
segment B overlaps with segment A. The total world-sheet area of the overlap is yl_B/2 and the average time a 
given infinitesimal segment (which carries on the information about B) stays on A is A, i.e. C2 — B /2 or similarly 
C4 = A/2. Thus, 



A = -HAB - HBA 



l-2v^ 



and 



The resulting operator 



H = 2v^ - SA ■ 



B = ^HBA - HAB 



l-2v2 5 d 



l-2v^ 



2 dA dB 



AB 



i-2v2 a d 



2 dB dA 



(A7) 



(A8) 



(A9) 



describes the evolution of correlation vectors in the transport equation Q under assumption of negligible back- 
reaction. Although we have only calculated the gravitational effects for Friedmann backgrounds it should be relatively 
easy to generalize the formalism to more complicated backgrounds, but is not immediately clear how to include back- 
reaction which lies beyond our semi-classical analysis. 



Appendix B: Longitudinal collisions 

Consider a longitudinal collision between two left-moving correlation vectors Bi and B3 and three right-moving 
correlation vectors Ai, A2 and A3 (see Fig. [3|. Under the string chaos assumption the directions of these correlation 
vectors are completely random but not all of such regions would emit simple loops. To distinguish the configurations 
that could lead to production of loops from those that cannot and to understand the aftermath of such coUis ons it 
is useful to plot a projection of unit vectors Ai, A2, A3, Bi, B3 on a plane orthogonal to B3 — Bi (see Fig. 4|. It 
is easy to check that for sufficiently large Ai, A3, Bi, B3 a loop of size ^ A2 would form whenever the 61,63 point 
lies inside of the triangle with vertices in Ai, A2 and A3. And if so, then the vectors Ai and A3 must be correlated 
with each other. (Such correlations are represented on Fig. [3] with a finer line along the kink connecting AiQ!i and 
A3q;3.) Thus, in order for the string chaos assumption to remain valid we will approximate such loops with only four 
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FIG. 3: Longitudinal collision of two left-moving correlation vectors with three right-moving correlation vectors. 



correlation vectors Bi/3i, B3/33, A2 and Aiai + A-^a^ which could be described in terms of the distribution function 
^- . . 

It is convenient to analyze the transitions of infinitesimal segments though boundaries 
KL,LM,MN,NO,OP,PQ,QR,RS and ST (see Fig. [s]) with separate integral expressions. In particular the 
transitions across MN, NO, OP, PQ, QR would describe a flow of energy from f to f and transitions across 
KL,LM,RS,ST describe a flow of energy from / to /. The corresponding creation of infinitesimal segments along 
boundaries AdN, NO, OP, PQ, QR is given by the following integral: 

/ Yi da^dA.dMBA f{Ai,Bi)f{A2,B2)f{A3,B3) S(^^^^_^ ^ja,A,_p^B,))S(^^^^^ (Bl) 

yi=i,2,3 J 

'^(Bi-B2)'^(a2-l)'^(fe) , Sr X (ft^r^' 
————^) ClCA-aiAi-aaAa) 2^ O(B-ftBi) ft + " 



A, 

ij^] + 0(A-A2)0(B-/JiBi)-g^ 



where 



A1A2 J B1B2 



For example the term with ^(A-aiAi-asAa) describes transitions across MN, NO, PQ,QR boundaries and the term 
with (^(A-Aa) describes the transitions across OP boundary. Similarly, the corresponding creation of infinitesimal 
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FIG. 5: The boundaries K L, LM, MN, NO,OP, PQ,QR, RS and ST describe possible transitions of infinitesimal segments 
during longitudinal collisions. 
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segments on boundaries ifL, LM, RS, ST is given by 



yi=l,2,3 J 

/ ^(Bi-B,)^(a,^i)'^(fe) , sr x (^ R \^^\\ 

and the overall destruction of infinitesimal segments on all boundaries {KL^ LM, MN, NO, OP, PQ, QR, RS, ST) is 
given by 

n ^^"^^^Aid^idBJ /(Ai,Bi)/(A2,B2)/(A3,B3) (5(52,^,,3(a.A.-/3.BO)'5(E.=i,2.3(".^.-ftBO (B4) 

^ ^,A.A3ivf |^E/(a-a.^(b-b4i 



'5(A-A2)'5(B-B2) 



B2 



The three integrals (Bl B3 B4) describe longitudinal collisions of two left-moving correlation vectors with three 



right-moving correlation vectors, but similar expressions would also hold for collisions of two right-moving vectors with 
three left-moving vectors with all of A's interchanged with all of B's. If desired one can also include more complicated 
longitudinal collisions involving more than five correlation vectors which would correspond to higher order effects. 



